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Center Manifold Theory 
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where the eigenvalues of A; have zero real parts and the eigenvalues 
of Az have negative real parts. 


Rewrite x = f(x) in the new coordinates: 


y = Ay + 81(y,2) 
Z = Az + g0(y,Z) 


g;(0,0) = 0, (0,0) = 0, (0,0) =0,i =1,2. 


Theorem 1: There exists an invariant manifold z = h(y) defined in a 
neighborhood of the origin such that 


Reduced System: y = Ayy+ei(y,h(y)) ye R* 


Theorem 2: If y = 0 is asymptotically stable (resp., unstable) for the 
reduced system, then x = 0 is asymptotically stable (resp., unstable) 
for the full system x = f(x). 


* Licensed under a Creative Commons 
Attribution-NonCommercial-ShareAlike 
4.0 International License. 


Khalil (Section 8.1), Sastry (Section 
7.6.1) 


EE C222/ME C237 - SPRING’18 - LECTURE 5 NOTES 2 


Characterizing the Center Manifold 
Define w = z —h(y) and note that it satisfies 
oh 
cae a 1 nae (Ay +81 (v,2)) 


The invariance of z = h(y) means that w = 0 implies w = 0. Thus, the 
expression above must vanish when we substitute z = h(y): 


Aak(y) + aly, My) ~ 5° (Ary +su(¥-h(y))) =0 


To find h(y) solve this differential equation for h as a function on y. 


If the exact solution is unavailable, an approximation is possible. For 
scalar y, expand h(y) as 


h(y) = hoy? + +++ + hpy? + O(y?*!) 


where hy = hg = 0 because h(0) = $°(0) = 0. The notation O(y?*?) 
refers to the higher order terms of power p + 1 and above. 
Example: 

y= Ye 

2=>-2+0y* «£0 


This is of the form (2) with ¢1(y,z) = yz, g2(y,z) = ay*, Az = —1. 
Thus h(y) must satisfy 


Reduced System: y = y(ay* + O(y?)) = ay? + O(y*). 
If a < 0, the full systems is asymptotically stable. If a > 0 unstable. 
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Discrete-Time Models and a Chaos Example 


Cl 22) = Fixit) DE: aaj =f) BHU ss 
fe = f(x*) =x* (“fixed point”) 
Asymptotic stability criterion: Asymptotic stability criterion: 
a a 
RA;(A) <0 where A= a - |A;(A)| <1 where A= a -_ 
f'(x*) <0 for first order system |f’(x*)| <1 for first order system 


These criteria are inconclusive if the respective inequality is not strict, 
but for first order systems we can determine stability graphically: 


Cobweb Diagrams for First Order Discrete-Time Systems 


Example: x,41 = sin(x,) has unique fixed point at 0. Stability test 
above inconclusive since f’(0) = 1. However, the "cobweb" diagram 


below illustrates the convergence of iterations to 0: 


Xn 


q 


. "THEIL, 
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Detecting Cycles Analytically 


fp=9a fg=p = f(r) =p fF@)) =4 


For the existence of a period-2 cycle, the map f(f(-)) must have two 
fixed points in addition to the fixed points of f(-). 


Period-3 cycles: fixed points of f(f(f(-))). 


Chaos in a Discrete Time Logistic Growth Model 


Xe = 11 — 2a) mn (3) 


Range of interest!:0<x<1 (Xn, >1 > Xn41 <0) 


r/4 


We will study the range 0 < r < 4 so that f(x) = r(1—x)x maps [0,1] 
onto itself. 


x* =0 and 


Fixed points: x =r(1—x)x => Sidi Sd, 


r <1: x* =0 unique and stable fixed point 


0 1 


r > 1: x =0 unstable because f’(0) =r >1 


j=) 
ray 

| 

In 

ray 
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Note that a transcritical bifurcation occurred at r = 1, creating the 
new equilibrium 


x*=1- 


2) eR 


Evaluate its stability using f’(x*) =r(1—2x*) =2-r. 
r<3 => |f'(x*)| <1 (stable) 
r>3 = |f'(x*)| > 1 (unstable). 


At r = 3, a period-2 cycle is born: 


FF (x) 

r(1— f(x)) F(x) 
=r(1l—r(1—x)x)r(1—x)x 
=rx(1—x)=r 


x 


=x =x) (l—r+re— 1") = x 
Factor out x and (x — 1+ 4), find the roots of the quotient: 


641 eVr-3)\(r +0) 


Pr or 
yx 
fee) Ff 
Be 
0 p I=eq 


This period-2 cycle is stable when r < 1+ V6 = 3.4494: 
E (F(x))) =f @)F() =f (PF @) = 4427-7 
P 


— 


J4+2r—r| <1 > 3<r<14+v6=3.4494 
At r = 3.4494, a period-4 cycle is born! 


“period doubling bifurcations” 
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m1 =3 period-2 cycle born 
12 = 3.4494  period-4 cycle born 
r3 =3.544 — period-8 cycle born 
14 = 3.564 — period-16 cycle born 


Too = 3.5699 


After r > fo, chaotic behavior for a window of r, followed by win- 
dows of periodic behavior (e.g., period-3 cycle around r = 3.83). 


Below is the cobweb diagram for r = 3.9 which is in the chaotic 
regime: 


